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Abstract 

This  report*  is  a  comparative  study  of  the  GPS-based  parameter  estimates  of  the  precise 
local  clock  linear  time  model.  We  estimate  the  time  error  x  and  fractional  frequency  offset  y 
employing  the  finite  impulse  response  (FIR)  digital  filters,  namely:  the  simple  moving  average 
(MA),  the  second  order  low  pass  (LP),  and  the  optimally  unbiased  (OU).  The  estimates  are  pro¬ 
vided  analytically,  and  we  show  that,  in  terms  of  a  minimum  root-mean-square  error  (RMSE), 
the  simple  MA  is  best  to  use  for  x  estimates  when  0  <  y  <  yx ,  the  second-order  LP  FIR  is  best  if 
yt  <  y  <  y2  >  alul  the  OU  is  best  for  y2  <  y ,  where  y  and  y2  are  determined  constants.  Fi¬ 
nally,  we  present  an  optimal  FIR  algorithm  for  the  precise  clock  linear  time  error  model  and 
show  that  of  all  the  estimators,  including  the  Kalman,  the  presented  algorithm  yields  minimum 
RMSE  of  the  fractional  frequency  offset. 


INTRODUCTION 

GPS-based  measurement  of  the  local  clock  performance,  the  time  error  and  the  frequency  offset, 
is  now  a  key  tool  in  timekeeping  [1].  Herewith,  it  is  shown  that  fast  and  accurate  performance  estimating 
is  constrained  by  a  large  variance  in  the  induced  GPS  noise.  The  standard  deviation  of  the  GPS  noise  us¬ 
ing  commercially  available  GPS  time  receivers  is  about  30  ns,  can  reach  10-20  ns  [2],  and  may  be  im¬ 
proved  by  removal  of  systematic  errors  to  3-5  ns  [2,3].  Even  so,  the  time  transfer  random  error  has  a 
much  higher  variance  than  that  of  a  precise  clock  and  the  measured  data  cannot  be  used  straightforwardly 
to  steer  the  clock  error.  Certainly,  this  is  a  stochastic  estimation  problem,  which  is  solved  in  practice  with 
different  filtering  techniques  requiring  processing  of  the  measurement  for  hours  and  even  days.  With  this, 
traditionally,  the  time  error  model  of  a  local  clock  is  associated  with  the  finite  polynomial  [4],  which  in 
discrete  time  is 


D  .22~ 
x„=x0+y0^n  +—An  +xn 


(1) 


where  x0  is  an  initial  time  error,  y0  is  an  initial  fractional  frequency  offset,  D  is  a  linear  fractional  fre¬ 
quency  drift  rate  (basically  representing  oscillator  aging),  x(t )  is  a  random  time  error  component, 
A  =  tn-  tn_x  is  a  sample  time,  t„  is  discrete  time,  and  n  =  0,  1 ,  ... 
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Timekeeping  employs  several  types  of  precise  clocks,  namely  atomic  (hydrogen,  cesium,  and  rubidium) 
and  crystal.  The  difference,  respecting  the  model  (1),  is,  first,  in  the  magnitudes  of  y0  and  D,  which  are 

strongly  limited  by  physical  nature  and,  second,  in  sensitivity  to  environment.  In  atomic  clocks  the  aging 
rate  D  normally  is  small,  the  proper  term  in  (1)  is  negligible  for  several  days  or  even  weeks,  and  the  time 
error  measurement  behaves  like  a  linear  noisy  process.  Respecting  this  property,  the  clock  performance 
estimators  in  modem  timekeeping  systems  are  usually  designed  based  on  either  the  two-state  Kalman  fil¬ 
ter  [5]  or  the  FIR  filter  with  a  constant-weight  ( simple )  MA  or  averaging  [3]. 

Of  all  known  filters,  the  simple  MA  is  best  in  terms  of  minimum  produced  noise  [6].  The  estimate  bias 
(inaccuracy),  however,  is  great  here.  Aiming  to  improve  accuracy,  we  have  recently  designed  a  new  MA 
filter  [7]  with  optimally  compensated  bias.  Because  the  simple  MA  is  optimal  in  terms  of  minimum  noise 
[6]  and  the  new  filter  [7]  is  optimal  in  terms  of  minimum  bias,  all  the  other  possible  FIR  filters  are  inter¬ 
mediate.  Among  them,  the  second-order  LP  filter  is  the  proximate  simple  rival. 

In  this  report  we  focus  attention  on  the  aforementioned  FIR  filters,  namely:  simple  MA,  second-order  LP, 
and  OU,  aiming  to  compare  their  performance  (estimates  of  time  error  and  fractional  frequency  offset)  for 
the  linear  time  error  model,  D  =  0,  perturbed  by  GPS  noise.  To  determine  the  tradeoff,  we  first  analyti¬ 
cally  provide  the  estimates  and  RMSE  for  each  filter.  Then  we  compare  the  filter  performance  for  the 
typical  magnitudes  of  the  fractional  frequency  offset  y0  and  GPS  noise  related  to  the  precise  clock.  Fi¬ 
nally,  we  present  the  optimal  FIR  estimation  algorithm  for  the  linear  time  error  model. 


TIME  ERROR  ESTIMATING  WITH  FIR  FILTERS 


In  the  time  error  measuring  system,  the  GPS  receiver  induces  into  the  signal  (1)  an  extra  noise  vn ,  which 
is  Gaussian  [8]  with  mean  zero  and  of  known  constant  variance  .  The  GPS-based  measurement  is 
combined  as  an  additive  sum 

z  —  x  +  v  .  (2) 

n  n  n  v  ' 

So  long  as  the  noise  variance  Gv  is  much  greater  than  that  of  the  model  (1),  the  term  xn  in  (1)  may  be 
neglected.  Then  for  the  linear  case,  D  =  0,  the  time  error  function  (1)  simplifies  to 

Xn  =  *0  +  ToAn  ’  (3) 


in  which  aging  may  be  accounted  for  by  updating  x0  and  y0  for  the  start  point  of  the  estimation  interval. 
In  view  of  (2),  the  FIR  estimate  of  a  time  error  is 


+  v(,W, 


(4) 


where  Wt  is  a  weighting  function  (kernel  or  impulse  response);  z N(ri)  =  [znzn_l...zn_N+ is  the  measure¬ 
ment  data  vector,  of  dimensions  TVxl ;  and  Wv  =  [IF0(A),fF1(A),...,fFjV_1(A)]r  is  the  filter  weight  matrix, 
of  dimensions  Ax  1 .  In  accordance  with  (2),  the  measurement  may  be  written  as  z N(n)  =  xN(n)  +  v  v  ( n ) , 
where  xN{n)  =  [xnxn_l...xn_N+l'\T  is  the  time  error  data  vector,  of  dimensions  Ax  1,  and 
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vN{ri)  =  {vnvn_x...vn_N+x\T  is  the  GPS  noise  data  vector,  of  dimensions  Vx  1 .  It  then  may  readily  be  shown 
that  for  the  linear  trend  (3),  the  estimate  (4)  becomes 

*n  =  xo  +  VoA(«  -  kI- )+  wr„ ,  (5) 


where  =  [0  1  ...N  -l]r  is  a  vector  of  integers  of  dimensions  /Vx  1 .  The  noise  wxn  =  v 'v  Wv  is  formed 
by  a  weighted  integration  of  the  white  Gaussian  origin  vn ,  presenting  samples  of  the  weighted  Wiener 
process.  Let  us  note  that,  after  integration,  the  low  frequency  components  of  the  noise  wxn  inherently 
dominate,  promising  some  estimation  trouble  caused  by  temperature. 

The  first  derivative  of  (5)  produces  an  estimate  of  the  fractional  frequency  offset;  this  is 


y,  =  ^(w  =  Vo  +  'v 


(6) 


in  which  increments  of  samples  of  the  weighted  Wiener  process 


1 

=  tK 


■  w 


-J 


x(n-l) 


=  —  [v T  (n)  -  vTN(n  -  1)]WV  rcPrcscnt  the  white  Gaussian  noise  (Appendix  A).  To  evaluate  the  noise  val- 
A 

ues  of  v n ,  wxn  and  vm  just  by  simple  MA,  Figure  1  shows  an  example  for  av  =30  ns,  N  =30,  and 
A  =  lsec.  It  then  insures  that  the  white  Gaussian  noise  v  of  the  frequency  estimate  (6)  is  appreciably 
lower  than  the  noise  wxn  of  the  time  error  estimate  (5). 


Let  us  now  evaluate  the  estimate  error  as  follows: 


£  „  =x„ 


(7) 


The  RMSE  is  then  calculated  by  Evn  =  ^E[e2xn ]  =  ^/Ar^  +  ~o\xn  ,  where  Axn  =E[&n~\  is  a  bias  (inaccu¬ 
racy),  and  a\xn  =  E[{ exn  -  Axn  )2  ]  is  a  variance,  representing  precision  of  the  filter.  Use  (3),  (4),  and  (7), 
take  into  consideration  that  x^W N  =  W2 x  v ,  and  write 


El„=E[(y0AkTNWN-vTNWN)2] 

=  y^A\klWN)2+E[(yTN\VN)2]-2y0AkTNWN\VTNE[yN].  (8) 

Because  the  noise  vn  is  mean  zero,  E  [ v  v  ]  =  0  ,  and  its  variance  is  constant,  then,  first,  the  last  term  in  (8) 
tends  to  zero,  and,  second,  the  transformation  i?[(v^  WN  )2  ]  =  O  2  W  v  W  2  brings  (8)  to 

=  ToA2(k^WjV)2  +a2W^Ww.  (9) 

Now  note  that  the  first  term  in  (9)  is  a  constant  bias  of  the  second  order,  which  means  that 

M  =  t0  Ak^W^,  (10) 
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and  the  second  term  in  (9)  is  in  turn  the  noise  constant  variance  of  the  filter: 


a2  =o2WrW 

u  ex  u  v  r y  N  TV  ' 


(11) 


Following  the  same  scheme  provides  the  error  estimate  for  the  fractional  frequency  offset  (6);  this  is 
E2„  =£[(f0  -j)J2],  and  write 


It  follows  from  (12)  that  first  the  estimate  of  the  fractional  frequency  offset  is  unbiased,  A yn  =  0,  and, 
second,  the  standard  deviation  of  the  estimated  noise  is  reduced  by  A  times,  being  dependent  on  the  filter 
weight.  It  in  turn  means  that  the  RMSE  (12)  is  in  fact  a  variance  of  the  estimate  error,  E2  =  a:,  .  Now 
let  us  examine  the  selected  FIR  filters  for  the  estimate  errors. 


SIMPLE  MA 


A  simple  MA  estimate  is  obtained  via  (5)  with  a  constant  weight  W  v  =  N~ll N ,  where  I  w  =  [1  1  ...if 
is  a  unit  matrix  of  dimensions  Ax  1,  and  with  a  delay  on  the  transient 

0  =  A(A-1).  (13) 

Substituting  the  above-mentioned  constant  weight  into  (5)  yields  the  time  error  estimate 


x„  =x0+y0A 


N- 1 


(14) 


where  w  is  a  noise  of  the  estimate  free  comments  for  (5)1  with  a  variance  a2  =  a2  /  N .  It  follows 

xn  L  \  /  j  vitx:  v 

from  (14)  that  the  bias  (10)  is  Ax  =  O.5>’o0 ,  and  the  RMSE  becomes 


E2=O.25y202  +a;i-.  (15) 

N 

The  commonly  known  conclusion  follows  from  (15):  the  estimate  variance  is  reduced  by  N  times  and  the 
bias  is  50%  for  the  linear  case.  The  simple  MA  estimate  of  a  fractional  frequency  offset  is  given  by  (6), 
where  the  estimate  noise  v  has  a  variance  a2  ,  which  in  view  of  (12)  is  the  second  order  of  the  RMSE 
of  the  estimate  (6)  (see  Appendix  A.  A): 


A2  N2  ' 


(16) 
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SECOND-ORDER  LP  FILTER 

The  second-order  LP  filter  is  superior  to  the  simple  MA  in  terms  of  bias.  Its  truncated  impulse  response 
is  calculated  by 


W„  = 


k-e 

x 

0 


0  <  n  <  M  - 1 
otherwise 


(17) 


where  X  is  a  time  constant;  M  represents  a  reasonable  weight  length,  which  is  normally  taken  to  be  3x ; 


and  k  = 


N-l 

t2> 


is  an  adjusting  coefficient  tending  W„  to  unit  area.  To  compare  errors,  it  is  neces¬ 


sary  to  tune  the  filter  for  the  transient  (13)  of  a  simple  MA,  which  is  X  =  A(N  —  1)  /  3 ,  and  set  M  =  N. 
Accordingly,  the  level  y  at  the  truncated  point  of  the  impulse  response  is  calculated  by 


y  =  e 


— (JV- 1) 


=  e~3  =  0.05  . 


(18) 


Substituting  (17)  into  (5)  yields  the  time  error  estimate.  The  bias,  after  routine  manipulations  using  of  the 
arithmetico-geometric  progression  and  q  =  e~AA  T- 1 ,  becomes 


k  A2To  (N  - 1)  L  g-A/T[i-g-A0V-|)/Ti 
xd-e-A/T)|  (A-D(i-e-A/T) ; 


which  taking  into  account  (13)  for  an  arbitrary  y  reduces  to 

A£  =  Toe[1_Z)(e’Y)]> 


(19) 


(20) 


where  an  auxiliary  function  Z)(0,  y)  becomes  constant  for  the  great  number  of  samples  (observe  that 
GPS-based  measurement  involves  basically  more  then  100  points  to  achieve  the  accurate-enough  esti¬ 
mate)  and  for  y  =  0.05  yields 


0(0,  y) 


Alny 

0(yA/e  -If 


0  +  A(1  -  y) 

0 


1«JV 

y=0.05 


0.684. 


(21) 


It  then  turns  out  that  (21)  simplifies  (20)  to  a  relation  Ar  =  O.316yo0  with  an  error  of  10%,  5%,  and  <1% 

for  9  <  /V,  1 5  <  /V,  and  1 00</V,  respectively.  That  means  that,  for  the  sample  time  of  A  =  100  s  and  say  5 
hours  of  averaging,  one  deals  with  N  =  180  samples,  so,  practically,  the  above-given  approximate  formula 
is  adequately  accurate. 


A  trivial  statistical  transformation  of  the  random  part  in  (5)  for  the  weight  (17)  produces  the  estimate 
noise  variance 


a 


2 

wx 


A2  0  £4  -2  — (JV-1-/J 

^  j=0 


(22) 
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-2— (JV-l-i)  -2 -i 

To  bring  (22)  to  the  final  form,  let  us  note  that  \e  T  =  \e  x  ar|d,  once  the  transient  is  over 

i= 0  1=0 

"  l 

with  N  <  i,  the  unlimited  sum  V  a = -  may  be  considered  instead.  It  brings  (22)  to  the  formula 

k= 0  1  —  <7 


_2  _2 

o . -  a, 


1 


wx  v  2 


x  l-e 


-2A/X 


which,  taking  into  account  (13),  produces 


=o;^G2(0,y) 


(23) 


in  which  6(0,7) ,  like  the  case  of  (21),  becomes  constant  for  1  «  N.  This  for  7  =  0.05  produces 


G2(9,7)  = 


Ain2  7 

0(1  —  72A/e )  i«iv 

V  *  7  7=0.05 

_2  A 


=  1.5.  It  then  may  be  shown  that  the  variance  (23)  is  calculated  by 


G  =1  .5g  —  with  almost  the  same  errors  as  in  the  case  of  (21). 

WX  V  0  v  2 


Based  upon  (5)  and  (20),  we  now  may  write  the  time  error  estimate  provided  with  the  second-order  LP 
filter,  which  is  xn  =  x0  +  t0A{/7  -  k(N  -  l)[l  -  D(0, 7)]}+  wxn ,  that  for  1  «  N  ,  7  =  0.05  ,  and  k  =  1 


simplifies  to 


xn  =x0+  ToA 


N-_ T 
3.16  y 


wr 


(24) 


where  the  variance  of  the  estimate  noise  wxn  is  given  by  (23). 

Respectively,  RMSE  of  the  estimate  (24)  with  account  of  (20)  and  (23)  is  obtained  with 

E2  =  yl Q2k2  [1  -  D(0, 7)]2  +  a2  —  G2  (0, 7)  and  simplifies  to,  as  in  the  case  of  (24), 

0 


E;=O.ly202  +G 


1.5 

AM 


(25) 


The  estimate  of  the  frequency  offset  has  the  common  form  (6),  where  the  noise  variance,  in  view  of  (12) 
and  Appendix  A.B,  is  equal  to  the  second  order  of  the  RMSE  of  the  estimate,  that  is: 


Cy  9 
A2  (fV-1)2  ' 


(26) 


OPTIMALLY  UNBIASED  MA 

The  OU  FIR  filter  was  designed  in  [7]  especially  for  the  linear  case  (3)  with  the  weighting  function 
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2(2 N  - 1)  -  6 i 
N(N  + 1)  ’ 

0, 


0<i<N-\ 


otherwise 


(27) 


Substituting  (27)  into  (5)  shows  that  bias  becomes  identically  zero,  and  then  the  time  error  estimate  is  de¬ 
termined  by 


=X0+  y 0Atl  +  Wxn  > 


where  the  noise  variance  is  equal  to  the  RMSE  of  the  estimate,  that  is: 

2  2(2# -1) 


E 1  =  a2  =c 

v  wx 


N(N  + 1) 


(28) 


(29) 


It  follows  from  (29)  that,  for  a  large  N,  the  standard  deviation  of  the  estimate  is  2  times  greater  than  that 
of  the  simple  MA.  Finally,  the  estimate  of  the  frequency  offset  (6)  is  accompanied  here  with  the  noise 
variance  and  RMSE  obtained  in  Appendix  A.C;  these  are  equal  to: 


2=  2  _  q;  4(5  A2  +4A  +  8) 
*  ^  “  A2  A2(A  +  1)2 


(30) 


A  COMPARATIVE  ANALYSIS  OF  THE  FILTERING  ERRORS 

Now  examine  all  three  FIR  filters  for  the  RMSE  of  the  time  error  and  fractional  frequency  offset  esti¬ 
mates,  relating  results  to  the  GPS-based  measurement  of  the  precise  clock  time  error. 

Estimating  the  Clock  Time  Error 

Figure  2  qualitatively  exhibits  trends  (15),  (25),  and  (29)  of  the  RMSE  for  the  variable  v0  and  with  the 
other  parameters  constant.  The  minimum  RMSE  of  each  filter  separates  values  of  y0  into  three  ranges. 
If  0  <  y0  <  y1,  the  simple  MA  yields  a  minimum  RMSE;  the  second-order  LP  filter  is  best  in  the  range 
yx  <  v0  <  y2  ;  and  the  optimally  unbiased  is  most  accurate  if  y,  <  y0  with  its  constant  error.  Solution  of 
an  equality  of  (15)  and  (24)  produces  the  first  cross-point  coordinate 


y  i  = 


a. 


-G2(0,y)-  — 
0  N 


0  ]j  O.25-[1-Z)(0,y)]2 


1  «N 
7=0.05 


(31) 


Making  the  same  calculations  for  (24)  and  (29),  we  obtain  the  second  point 


T2  = 


o, 


2(2N  - 1)  A 


0[1-O(0,Y)]  \  N(N  +  l)  0 


--G-(0,y)|  =  5.00ovl|— 


1  «N 
'  y=0.05 


(32) 
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By  definition,  the  minimum  possible  filtering  error  of  the  stationary  process,  y0  =  0 ,  is  provided  with 
simple  MA,  that  is: 


E 


xmin 


1 


(33) 


Substituting  (31)  into  (15)  or  (25)  either  produces  the  first  intermediate  RMSE  (Figure  1),  namely 


E  ,  =  1.355a  .  1  =  1.355a, J- 

Vtvm  v  e 


(34) 


Finally,  calculating  (25)  or  (29)  either  for  (32)  yields  the  second  intermediate  error: 


Ex2  ~  °v  Vam  2°v\e: 


(35) 


Consider,  for  example,  the  typical  case  of  the  GPS-based  measurement  of  a  clock  time  error;  this  is 
Cv  =  30  ns,  A  =  100  sec,  N=  865,  and  0  =  24  hours.  The  minimum  RMSE  (33)  calculates  for  y0  =  0 

to  be  Ermin  =  1  ns,  and  then  RMSE  rises  up  to  Evl  =  1.36  ns  (34)  with  v0  =  2.16  ■  10~14  ,  proving  that  in 
the  range  of  0  <  v0  <  2.16  ■  1014  the  simple  MA  is  best.  With  y0  =  5.9 1  ■  10  14  the  RMSE  reaches 
Ev2  =2  ns  (35),  and  then  the  second-order  LP  filter  is  best  in  the  range  of 
2.16-10  14  <3^0  <5.9110  14.  Assuming  5.91  ■  10~14  <  y0,  we  conclude  that  the  RMSE  remains 
E  x2  =  2  ns  using  the  optimally  unbiased  filter. 


ESTIMATING  THE  CLOCK  FRACTIONAL  FREQUENCY  OFFSET 

It  follows  from  (16),  (26),  and  (30)  that  the  RMSE  of  the  frequency  offset  estimate  is  y0  invariant,  and  its 

variance  is  appreciably  reduced  by  A2  and  A.  Furthermore,  as  the  simple  MA  produces  minimum  noise 
among  all  the  known  filters  [6],  including  the  Kalman,  it  also  provides  the  best  estimate  of  the  frequency 
offset  in  the  linear  case.  For  example,  the  above-considered  case  (Section  IV.A)  calculates  the  standard 
deviation  Ev  =  4.9-10~13  with  the  simple  MA,  Er  =  1 .04  - 1 0  12  with  the  second-order  LP  FIR,  and 

E  =  1.55-10  12  with  the  optimally  unbiased  FIR.  Yet,  in  the  reported  GPS-based  timekeeping  system 
[3],  the  measurement  noise  is  achieved  with  <yv  of  about  5  ns,  allowing  the  clock  frequency  steering  of 
10'13  or  better  for  an  averaging  time  of  0  =  lday.  The  consistent  estimate  is  calculated  via  (16)  as 
E,.  =  8.2  •  1(T14 . 


OPTIMAL  FIR  ALGORITHM 


It  now  follows  from  the  above-given  analysis  that  the  designed  FIR  filtering  algorithm  is  optimal  in  the 
sense  of  minimum  RMSE  for  the  precise  clock  online  performance  estimate;  this  is 


(36) 
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k  =  ^  fc(«) -**(”- ok* »  (3?) 

in  which  Ww  is  given  by  Ww  =  A-1 IN,  (17),  and  (27)  for  y0  <  yx,  y0  <  y2,  and  y2  <  y0,  respectively. 

Herewith,  the  weight  (17)  may  be  considered  for  the  other  LP  filter.  The  special  feature  of  the  algorithm 
is  that  it  inherently  produces  the  minimum  possible  RMSE  of  the  linear  model  frequency  offset  estimate 
(37)  among  all  the  filters,  including  the  Kalman. 


DISCUSSION 

GPS-based  measurement  of  the  time  error  of  a  precise  local  clock  looks  like  a  linear  noisy  function,  as  the 
clock  frequency  drift  rate  is  negligible  for  a  few  days  or  even  decades.  Respecting  this  property,  we  have 
set  D  =  0  in  (1)  and  examined  the  errors  of  the  three  FIR  filters,  which  seem  to  be  most  accurate:  the 
simple  MA  is  optimal  in  a  sense  of  minimum  produced  noise  [6],  the  optimally  unbiased  FIR  filter  is  op¬ 
timal  in  a  sense  of  a  minimum  bias  [71,  and  the  second-order  LP  FIR  filter  is  the  proximate  simple  rival 
with  its  intermediate  performance.  Depending  on  y0 ,  the  following  FIR  filter  produces  a  minimum 
RMSE  of  the  time  error  estimate  (Fig.  2),  namely:  the  simple  MA  is  best  in  the  range  0  <  y0  <  yt ;  the 
second-order  LP  is  best  in  the  range  yx  <  y0  <  y2 ;  and  the  optimally  unbiased  yields  a  constant  best  es¬ 
timate  if  y2  <  v0 .  The  estimation  filtering  algorithm  (36)  and  (37)  is  optimal  for  the  linear  time  error 
noisy  function  in  a  sense  of  minimum  RMSE. 


APPENDIX  A:  ESTIMATE  VARIANCE  OF  THE  FRACTIONAL  FRE¬ 
QUENCY  OFFSET 

Simple  MA  Estimating 


Consider  the  noisy  part  in  (6)  for  the  constant  weight  1/A;  this  is 

1 

AN 


A,,  =^-[v«00-v^0'-l)]w 


N- 1 


Z(r«  ,  i) 


=  TT7(k,  -  V,v)> 
AN 


(Al) 


in  which  increments  vn  —  vn_N  of  the  white  Gaussian  noises  form  the  white  Gaussian  noise  (Al)  of  the 
variance 


2  _  1 

^  ~  A2A2 


[(a  -  Vn-N  )2  ]  =  TT77T  [£(V«  )  +  E(Vn-N  )-  2E(NnVn-N  )]• 


A  A 


Observe  that  samples  of  the  shifted  stationary  white  Gaussian  noise  are  uncorrelated,  E{ynvn_N)  =  0  ,  and 
al  =  E( \>l )  =  E(v2n_N ) .  Then  pass  to  the  formula  (16). 
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Second-order  LP  FIR  Estimating 

Consider  the  noisy  part  in  (6)  for  the  weight  ( 1 7)  and  T  =  A(N  —  1)  /  3 ,  this  is 


v  = 

yn 


U  N- 1 


3 

'am' 


A(A^-l)~o 


3& 


A(A-l) 


r  3 

X 

3 

r  3 

X 

v  —v  , 

n  n- 1 

\-e^ 

L 

-v  e  ^ 

Vn- 2^ 

\-eN-' 

V 

J 

V 

7 

_ 

3 

Observe  that  for N»  l  the  decomposition  holds  true  of  e  N~l  =  1  —  3 /(N  —  1) .  Then  neglect  all  the  noise 
samples  reduced  by  N-l  in  the  brackets  and  write 


v  = 

yn 


3  k  ,  x  3  k 

-(v„+0.05v„_J  = - v„. 


A(A-l)  A(N  -V) 

The  noise  (A2)  is  white  Gaussian,  and  (A2)  readily  produces  the  variance  (26)  for  k  =  1 . 


(A2) 


Optimally  Unbiased  FIR  Estimating 

Consider  the  noisy  part  in  (6)  for  the  OU  weight  (27),  this  is 

2 


[  N- 1 


V  = 

yn 


AN(N  +  l) 

2  N~' 


X  [(2  A  - 1  -  3/)vv,  -  (2 N  - 1  -  3i ] 


AN(N  +  l)U 
2 


S  [(2iV  “  ^(v,-  -  V,-i )  -  3/(v„_,  -  v„_M )] 


AN(N  +  l) 


(2N  —  l)(vn  -  vn_N )  +  3 


JV-1 


r-„  +(A-l)v„_,v 


2 

AN 


3N 


„  N-2 

2v  H - V  ,, - 

"  A  +  l  A  +  l 


(  I  H  A 

— y  v  . 


Note  that  the  expression  in  parenthesis  of  (A3)  is  zero-mean  of  the  white  Gaussian  noise,  then 

2  f 


Vyn  AN 


2v„ 


N-2  r  ^ 
N  + 1  V"~r 


(A3) 


(A4) 


In  a  view  of  (Al)  the  noise  (A4)  is  white  Gaussian.  The  samples  vn  and  vn_N  are  uncorrelated,  so  the 
variance  of  the  noise  (A4)  is  readily  provided  in  a  form  of  (30). 
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Fractional  frequency  offset  _y0 


Figure  2.  RMSE  of  the  time  error  estimates  provided  with  the  FIR  filters  for  different  v0 


